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ABSTRACT 

Visibility  measurements  with  Michelson  interferometers,  particularly  the  measurement  of  fringe  contrast,  are 
affected  by  various  atmospheric  and  instrumental  effects,  all  of  which  reduce  the  measured  contrast.  To  compen¬ 
sate  for  this,  stars  with  known  or  predictable  diameters  (calibrators)  are  observed  so  that  the  overall  reduction 
in  the  visibility  can  be  measured.  Objects  with  the  smallest  possible  diameters  are  preferred  as  calibrators,  since 
the  predicted  visibilities  become  less  sensitive  to  any  uncertainties.  Therefore,  unreddened,  early  type  stars  are 
usually  chosen  if  they  are  available  because  they  are  relatively  bright  for  a  given  angular  diameter.  However, 
early  type  stars  bring  additional  complications.  Rapid  rotation,  common  with  these  stars  can  cause  variations 
in  the  visibility  amplitudes  due  to  oblateness  and  surface  brightness  asymmetries  that  are  larger  than  implied  by 
the  usual  error  estimates.  In  addition,  rotation  can  introduce  significant  phase  offsets.  Using  Roche  models,  von 
Zeipel  theory,  and  the  observed  constraints  of  V,  B-V,  and  vsin  i,  it  is  possible  to  put  limits  on  the  size  of  these 
effects  and  even  estimate  the  distribution  of  possible  visibilities.  To  make  this  easily  available  to  the  community, 
we  are  in  the  process  of  creating  a  catalog  of  possible  calibrators,  including  histograms  of  the  visibilities,  cal¬ 
culated  for  configurations  used  at  a  number  of  observatories.  We  show  the  examples  of  several  early  type  stars 
which  are  potential  calibrators  using  parameters  appropriate  for  the  Navy  Prototype  Optical  Interferometer. 

Keywords:  Calibrators,  Interferometry,  Early  type  stars 

1.  INTRODUCTION 

Michelson  interferometers  in  principle  measure  the  complex  visibility,  which  is  the  Fourier  transform  of  the 
object’s  surface  brightness  integrated  perpendicular  to  the  projected  direction  of  the  baseline1 .  However,  the 
usual  interferometric  observables  are  the  squared  visibility  amplitude,  an  unbiased  estimator2  and,  if  three  or 
more  apertures  are  available,  the  closure  phase,  which  is  free  of  atmospheric  effects. 

Various  instrumental  and  atmospheric  effects  act  to  reduce  the  measured  visibility  amplitude  and  can  add 
offsets  to  the  phases.  These  can  be  corrected  by  observations  of  calibrators3  .  The  squared  visbility  amplitude, 
V2,  of  a  science  target  is  then  obtained  as 


V2  = 


Meal 


K2al 


(1) 


where  /i  and  /ucai  are  the  measured  coherence  factors  of  the  target  and  the  calibrator  respectively,  and  the  V2al 
is  the  squared  visiblitity  amplitude  of  the  calibrator,  for  which  only  estimates  are  usually  available.  Phases  can 
be  similarly  adjusted,  usually  by  additive  corrections. 
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Stars  having  the  smallest  possible  diameters  are  generally  chosen  as  calibrators  since  the  predicted  visibilities 
are  near  maximum  and  hence  less  sensitive  to  uncertainties  in  the  assumed  diameters,  which  are  of  the  order  of 
a  few  percent  when  estimated  by  surface  brightness-angular  diameter  relations3,4.  However,  most  stars  hotter 
than  F5  rotate  rapidly  and  rotation  can  introduce  significant  oblateness  and  surface  brightness  asymmetry  in 
the  apparent  disks5  which  in  turn  can  cause  variations  in  the  inferred  diameters  in  addition  to  what  is  implied 
by  the  usual  error  estimates.  In  addition,  asymmetry  can  cause  a  significant  triple  phase  signal. 

Here,  we  investigate  how  rotation  might  affect  the  use  of  the  B-F  stars  as  calibrators.  We  estimate  the 
magnitude  of  the  rotation  effects  using  von  Zeipel  theory5  and  the  constraints  by  provided  by  measurements 
of  the  magnitude  and  color,  V  and  B-V,  the  projected  equatorial  velocity,  vsini,  and  the  parallax,  tt.  In  the 
next  section  we  describe  the  gravity  and  limb-darkened  Roche  spheroid  model  of  a  calibrator.  We  introduce  the 
instrumental  parameters  and  the  criteria  for  selecting  B-F  stars  as  calibrators  in  section  3.  Then  we  describe  the 
calculation  of  the  simulations  in  section  4.  Lastly,  we  describe  the  catalog  of  the  calibrators  we  will  produce. 

2.  A  MODEL  FOR  EARLY  TYPE  STARS  AS  CALIBRATORS 

The  theory  of  rotating  stars  was  studied  by  von  Zeipel5  in  1924.  He  showed  that  a  rotating  star  adopts  the 
figure  of  a  Roche  spheroid  whose  effective  temperature  is  proportional  to  the  fourth  power  of  the  effective  gravity 
(Teff  oc  (<7e//)J)  under  the  assumptions  of  radiative  equilibrium  and  solid  body  rotation.  This  theory  has  been 
tested  not  only  in  close  binary  systems  but  also  through  interferometric  observations  of  isolated  stars6-8  .  When 
calculating  observables  with  this  model,  it  is  assumed  that  plane-parallel  atmosphere  describes  the  atmosphere 
locally.  We  use  the  ATLAS9  models  parametrized  by  Van  Hamme10  .  Since  these  are  hot  stars  we  have  adopted 
the  square  root  parameterization  of  limb-darkening11  . 

3.  SELECTION  OF  CALIBRATORS 

3.1.  Instrumental  Parameters 

Whether  rotation  introduces  important  asymmetries  depends  critically  on  combination  of  the  wavelengths  and 
baselines  used.  We  focus  here  on  the  Navy  Prototype  Optical  Interferometer12  (NPOI)  which  has  baselines  up 
to  80  m.  We  adopt  500  nm  as  the  nominal  wavelength  for  calculating  the  squared  visibilities. 

3.2.  Selection  of  Calibrators 

Table  1.  Examples  of  potential  calibrators  :  risini  is  in  units  of  km/s.  Angular  diameter,  Qbem  and  parallax,  7r  are  in 
mas.  The  angular  diameter  was  calculated  using  the  Barnes-Evans-Moffet  relation4  . 


HR 

Name 

RA 

DEC 

SpType 

V 

B-V 

vsini 

Obem 

7 r 

1673 

68  Eri 

050843.6 

-042722 

F2V 

5.12 

0.44 

10 

0.606 

39.99 

3974 

21  LMi 

100725.8 

+351441 

A7V 

4.48 

0.18 

148 

0.546 

35.78 

4733 

14  Com 

122624.1 

+271606 

FOp 

4.95 

0.27 

227 

0.505 

11.92 

8351 

51/i  Cap 

215317.8 

-133306 

F1III 

5.08 

0.37 

87 

0.554 

36.15 

Our  input  list  of  potential  calibrators  satisfies  the  following  criteria:  (1)  luminosity  class  III,  IV,  and  V 
stars  in  the  color  range  B-V<  0.45  as  listed  in  the  Bright  Star  Catalogue13,  (2)  declinations  above  —15°,  (3) 
angular  size  ( Obem )  estimated  using  the  Barnes-Evans-Moffet4  relation  between  0.5  mas  and  1.5  mas,  and  (4) 
distance  closer  than  100  pc.  The  last  eliminates  significant  reddening.  We  then  pruned  the  list  of  any  stars  with 
companions  with  separations,  p  <  1”  and  magnitude  differences,  Am  <  5.0.  However,  unless  second  spectra 
were  reported,  we  retained  those  listed  or  suspected  as  spectroscopic  binaries.  Examples  of  potential  calibrators 
are  shown  in  Table  1 . 
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Table  2.  The  maximum,  minimum,  weighted  average,  ±<r  values  of  V2  of  the  sample  stars  with  NPOI 


HR 

Name 

Vmax 

Vmin 

V 2 

-  a 

+CT 

1673 

68  Eri 

0.656 

0.629 

0.630 

0.629 

0.630 

3974 

21  LMi 

0.688 

0.677 

0.683 

0.678 

0.687 

4733 

14  Com 

0.743 

0.607 

0.656 

0.644 

0.671 

8351 

51/lz  Cap 

0.701 

0.678 

0.683 

0.680 

0.686 

4.  SIMULATIONS  AND  CALCULATIONS 

4.1.  Procedure  of  Simulations 

Roche  models  require  several  parameters.  The  initial  inputs  are  inclination  angle,  i,  mass,  angular  size,  6  (or 
equivalently  the  radius,  R,  through  the  parallax,  7r),  effective  temperature,  Teff,  and  the  dimensionless  angular 
velocity,  to  (=  Q/Qb  where  tt  is  the  angular  velocity  and  Qb  is  the  angular  velocity  at  break-up).  We  obtain 
a  mass  by  finding  the  model  in  the  Padova  grid  which  reproduces  the  luminosity  and  polar  radius.14  During  a 
simulation,  the  model  is  required  to  satisfy  to  the  observed  V,  B-V  and  vsin  i.  That  is,  we  iterate  until  we  match 
the  projected  rotational  velocity  (i/sini)  by  varying  the  angular  velocity  parameter  (w)  or  inclination  angle  (i), 
or  both.  In  the  process  the  program  must  also  iterate  for  consistent  values  of  the  polar  values,  Tef/.p  and  9P,  to 
match  V  magnitude  and  B-V,  respectively.  The  models  then  allow  us  to  predict  visibilities  as  a  function  of  the 
assumed  Position  Angle  of  the  rotational  axis. 

In  the  process  of  these  simulations  we  have  found  that  calculated  polar  gravity,  log  gp,  increases  monotonically 
as  the  inclination  angle  decreases.  The  result  is  that  for  many  candidates  the  potential  range  of  the  inclinations 
is  terminated  well  before  reaching  break-up  angular  velocity  ( to  =  0.999)  because  the  predicted  gravities  would 
exceed  those  of  the  predicted  Zero  Age  Main  Sequence  (ZAMS). 

Squared  visibilities  can  be  calculated  at  any  inclination  angle  and  a  Position  Angle  allowed  by  the  mentioned 
constraints  of  the  observed  quantities.  The  allowed  inclination  angles  are  between  i  =  90°  and  the  ’’critical 
angle”  if  limited  by  ZAMS  constraints  on  the  gravity.  However,  if  there  are  no  restrictions  of  the  gravity,  we  do 
simulations  up  to  the  break-up  angle.  The  squared  visibilities  are  sampled  at  12  evenly  spaced  Position  Angles 
and  20  evenly  spaced  inclinations  in  the  allowed  range  of  each  star. 

4.2.  Statistical  Calculations 

Using  the  visibilities  calculated  over  the  range  of  inclinations  and  Position  Angles,  we  generate  the  histograms 
describing  the  probability  distribution,  and  calculate  the  maximum,  minimum,  average,  and  95%  confidence 
levels  of  V2.  The  only  assumption  here  is  that  the  rotational  axis  is  randomly  oriented15  over  the  range  of 
allowed  inclinations  and  Position  Angles. 

The  histograms  are  constructed  by  dividing  the  range  of  squared  visibilities  into  ten  bins.  The  probability, 
Pk  of  squared  visibilities  being  the  k  th  bin  is  then  given  by, 


Pk  =  P(V2  <v2<  V2+1)  =  ^!*sintl ,  (2) 

z2j  sinij 

10 

P  =  Y,Pk  =  l,  (3) 

fc= i 

where  N  and  are  the  total  number  of  inclination  and  position  angle  samples  and  the  number  that  satisfy 
V2  <V2<  V2+l  respectively.  And  ij  and  ii  are  the  inclination  angles  at  the  sample  points. 


Proc.  of  SPIE  Vol.  6268  626848-3 


HR1673 


Figure  1.  The  histogram  for  HR  1673:  This  star  is  a  slow  rotator  with  nsini  =  lOkm/s.  In  spite  of  its  slow  rotation  and 
a  narrow  range  of  ±cr  the  histogram  displays  an  extended  tail,  albeit  of  low  probability.  These  extreme  visibilities  result 
from  the  fact  that  star  could  be  seen  nearly  pole  on.  The  star  should  be  used  as  a  calibrator  with  some  caution. 


The  weighted  average  of  the  squared  visibilities,  V 2,  of  the  star,  corresponding  to  the  inclination  angles,  is 
obtained  from 


V2  = 


EN  X  7-0 

j  VjSinij 
N 


(4) 


where  V '2  is  the  squared  visibility  at  the  j  th  sample.  Finally,  we  tabulate  squared  visiblities  of  the  95  %  level  of 
significance,  indicated  as  ”±cr”,  as  well  as  the  maximum  (V^ax)  and  minimum  (V2nln)  values. 


5.  RESULTS 

We  limit  our  discussion  here  to  visibility  amplitudes.  Estimates  in  inferred  angular  diameters  by  the  surface 
brightness-angular  diameter  relations  are  uncertain  by  a  few  percent.  Adopting  the  3%  uncertainty  suggested  by 
Muzurkwich  et  al.3  in  a  star  with  a  unform  disk  angular  diameter  of  0=0.55  mas,  measured  with  a  80  m  baseline 
at  550  nm  experiences  an  uncertainty  in  V2  =  0.63  of  about  0.02.  The  effects  of  rotation  on  the  visibilities  can 
be  substantially  larger  than  this. 

To  illustrate  the  range  of  results,  we  consider  4  cases  listed  in  Table  1.  They  include  a  slow  rotator  (HR 
1673),  an  intermediate  rotator  (HR  8351),  and  two  moderately  fast  rotators  (HR  3974,  HR  4733).  Note  that 
by  usual  standards  all  would  at  first  glance  be  good  calibrator  choices,  with  angular  diameters  nominally  in  the 
range  0.5-0. 6  mas. 

Table  2  presents  the  results  of  the  simulations  where  V£in,  V^ax,  V2,  and  ±cr  values  of  the  squared  visibilities 
of  each  star  are  given.  The  histograms  derived  from  these  calculations  are  shown  in  Fig  1  through  Fig  4  which 
show  the  probability  distribution  of  the  V 2  values  according  to  the  allowed  range  of  inclinations  and  Position 
Angles.  The  bin  size  for  the  histograms  is  a  tenth  of  the  range  V^ax  to  V^in  if  this  is  larger  than  0.2  percent  of 
V2,  and  otherwise  fixed  at  0.2%  of  the  average  if  the  range  is  less  than  2%  of  V2. 

Stars  with  small  projected  rotational  velocities,  predicted  to  be  or  near  the  ZAMS  such  as  will  show  no 
significant  variation  in  their  visibilities.  These  are  ideal  calibrators  and  we  have  identified  a  number  in  our 
catalog.  On  the  other  hand,  the  slow  rotator,  HR  1673,  illustrated  in  Fig  1,  shows  an  extended,  low  probability 
tail  in  the  squared  visiblilty.  Even  though  it  is  likely  to  have  the  nominal  visibility,  there  is  a  finite  chance  of  a 
4  %  error.  In  the  simulations  for  this  star  the  inclination  was  cut  off  at  an  inclination  of  i  =  2.45°,  with  smaller 
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Figure  2.  The  histogram  for  HR  8351:  This  star  is  an  intermediate  rotator  with  rsini  =  87km/s.  The  histogram  shows 
the  rotational  effects  on  its  visibility.  This  star  can  be  a  calibrator,  but  the  variation  of  squared  visiblities  should  be 
properly  considered  because  of  the  extended  tail  of  possible  V2  values.  The  abscissa  covers  the  entire  range  of  visibilites. 


HR4733 


Figure  3.  The  histogram  for  HR  4733:  This  star  is  a  fast  rotator  with  rsini  =  227km/s.  As  expected,  the  rotation 
significantly  affects  the  visiblilites.  This  is  an  example  of  a  star  with  a  small  predicted  angular  diameter,  0.505  mas,  which 
could  nevertheless  produce  significant  calibration  errors. 
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Figure  4.  The  histogram  for  HR  3974  :  This  is  a  moderately  rotating  star  with  vsini  =  148km/s  which  is  on  or  very 
near  the  ZAMS.  This  can  be  a  good  calibrator  because  of  the  small  predicted  range  in  V2 .  However,  it  is  clear  that  its 
rotation  does  produce  some  uncertainty  in  its  visibility.  Here,  the  size  of  bin  is  the  0.2  percent  of  V2. 


inclinations  pushing  the  models  below  the  ZAMS.  Even  with  critical  rotation  excluded,  the  high  inclination  that 
are  accessible  can  generate  significant  uncertainties  in  the  visiblilities. 

Fig  2  shows  the  histogram  of  squared  visibilities  for  HR  8351,  which  has  an  intermediate  projected  rotational 
velocity  of  87km/s.  This  rotation  can  result  in  significant  range  of  visibilities  for  this  star.  Although  a  reasonable 
candidate  as  a  calibrator  as  with  HR  1631,  this  star  can  add  uncertainty  to  visibility  measurements. 

For  the  moderately  fast  rotator  HR  4733  we  show  how  its  rotation  affects  the  predicted  squared  visibilities  in 
Fig  3.  Significant  flattening  is  predicted  at  all  inclinations,  so  the  visibility  measurements  vary  significantly  over 
the  possible  combinations  of  inclinations  and  position  angles.  In  spite  of  a  small  predicted  angular  diameter, 
0.505  mas,  the  range  in  the  visiblilities  must  be  considered  calibrating  with  this  star. 

In  Fig  4,  we  show  the  visibility  range  for  HR  3974,  a  star  with  moderately  fast  projected  rotation.  HR  3974 
appear  to  lie  on  or  very  the  ZAMS  which  greatly  limits  the  predicted  range  of  visibilities,  allowing  it  to  be  a 
good  calibrator  in  spite  of  its  significant  projected  velocity.  The  double  horned  shape  of  the  resulting  histogram 
is  characteristic  of  stars  near  the  ZAMS  affected  by  significant  rotation. 

By  looking  through  the  variety  of  histograms,  one  can  identify  the  kinds  of  early  type  stars  that  can  be  good 
calibrators.  The  expected  cases,  early  type  stars  with  small  angular  diameter  and  small  projected  rotational 
velocities  are  not  guaranteed  to  be  good  calibrators.  However,  stars  that  can  be  identified  as  on  or  near  the 
ZAMS  are  generally  good  calibrators  because  of  their  narrow  visibility  distribution  even  with  significant  rotation. 

6.  CONCLUSIONS 

We  are  preparing  a  catalog  of  early  type  stars  (unreddened  and  hot  stars)  with  small  diameters  that  can  serve  as 
calibrators.  Generally,  for  these  stars  the  predicted  visibilities  are  less  sensitive  to  uncertainties  in  the  predicted 
diameters  and  they  are  relatively  bright  at  a  given  angular  size.  However  fast  rotation,  common  with  these  stars, 
causes  variations  in  the  inferred  visibilities  due  to  oblateness  which  might  be  larger  than  implied  by  the  errors 
estimated  for  the  angular  diameters.  Using  Roche  models,  von  Zeipel  theory,  and  the  observed  constraints  of  V, 
B-V,  and  usinz,  we  estimated  the  range  in  visibilities  that  might  be  induced  by  rotation.  Visibility  histograms 
of  some  representative  examples  illustrate  the  effects  of  rotation  and  proximity  to  the  ZAMS.  Early  type  stars 
with  small  diameter  not  affected  by  fast  rotation  are  not  guaranteed  to  be  good  calibrators.  We  have  also  found 
that  stars  on  or  near  ZAMS  can  also  be  reliable  calibrators,  even  with  moderate  projected  rotational  velocities. 
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